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Abstract 

We present an independent calculation of that part of the 0{a 2 s ) contribution to 
the decay width T(B — > X s ^y) which arises from the self-interference term of the elec- 
tromagnetic dipole operator Oj. Using a different method, we find complete agreement 
with a previous calculation. This NNLL contribution is an important ingredient for 
the complete NNLL prediction of T(B — > X s j) which will resolve the charm quark 
mass ambiguity appearing at NLL accuracy. 
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1 Introduction 



Rare B decays as flavor changing neutral currents (FCNCs) are loop-induced processes and 
highly sensitive probes for new degrees of freedom beyond the SM. But this indirect search 
for new physics signatures within flavor physics takes place today in complete darkness, 
given that we have no direct evidence of new particles beyond the SM. However, the day 
the existence of new degrees of freedom is established by the LHC, specifically when the 
mass scale of the new physics is fixed, the searches for anomalous phenomena in the flavor 
sector will become mandatory. In this context, the measurement of theoretically clean rare 
decays, even when found to be SM-like, will lead to important and valuable information 
on the structure of the new-physics models and to complementary information to the LHC 
collider data. 

Because physics effects beyond the SM in flavor observables seem to be rather small, 
as recent experimental data indicate, the focus on theoretically clean observables within 
the indirect search for new physics is crucial in order to detect specific patterns and to 
distinguish between various new physics scenarios. Inclusive rare B decays belong to this 
class of observables because they are dominated by perturbative contributions and non- 
perturbative effects are generally suppressed by AQ CD /m^ within the heavy mass expansion 
(for a review see [1]). 

The inclusive b — > 57 mode is still the most prominent rare decay. The stringent bounds 
obtained from this mode on various non-standard scenarios (see e.g. [2-9]) indicate the 
importance of this specific FCNC observable in discriminating new-physics models. It has 
been measured by several independent experiments [10-14], and the present experimental 
accuracy is already below the 10% level. In the HFAG world average of these measurements, 
a common shape function for the extrapolation to low energies was used; for a photon energy 
cut E y > 1.6 GeV the average is given by [15]: 

BR[B -> X sl ] = (3.55 ± 0.24 + ° f ± 0.03) x 10" 4 , (1.1) 

where the errors are combined statistical and systematic, the systematic ones due to shape 
function, and due to the contamination. In the near future, more precise data on this 
mode are expected from the (upgraded) B factories and from the planned Super- 5 factories 
(see [16-18]). 

Clearly, the theoretical uncertainty should be reduced accordingly. After completion of 
the QCD NLL [19-28]. and the electroweak two-loop [29-32] computations, it was generally 
believed that the theoretical uncertainty of the branching ratio is below 10%. However, as 
first pointed out in [33], there is an additional uncertainty in the NLL prediction, which 
is related to the renormalization scheme of the charm-quark mass. In a recent theoretical 
update of the NLL prediction of the branching ratio of B — > Xg'y, the uncertainty related 
to the definition of m c was taken into account by varying m c /mb in the conservative range 
0.18 < mc/mb < 0.31, which covers both the pole mass (with its numerical error) value and 
the running mass m c (/i c ) value with [i c e [m c , m&]; for E^ > 1.6 GeV one then finds [34] 

BR[B X sl ] = (3.52 ± 0.32\ mc/mb ± 0.02| CKM ± 0.24| param . ± 0.14| scale ) x 10" 4 . (1.2) 

The only way to resolve this scheme ambiguity in a satisfactory way is to perform a systematic 
NNLL calculation, which according to a recent investigation [35] is expected to reduce this 
uncertainty by a factor 2. 
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Parts of the three principal calculational steps leading to the NNLL result within the 
effective field theory approach are already done: (a) The full SM theory has to be matched 
with the effective theory at the scale \x = [Mw, where nw denotes a scale of order m w or m t . 
The Wilson coefficients Ci(/j,w) om Y pick up small QCD corrections, which can be calculated 
in fixed-order perturbation theory. In the NNLL program, the matching has to be worked 
out at order a 2 . The matching calculation to this precision is already finished, including the 
most difficult piece, the three-loop matching of the operators Or ,8 [36]. (b) The evolution 
of these Wilson coefficients from /x = down to \x = fib then has to be performed with the 
help of the renormalization group, where /ij is of the order of mj. As the matrix elements 
of the operators evaluated at the low scale ^ are free of large logarithms, the latter are 
contained in resummed form in the Wilson coefficients. For the NNLL calculation, this 
RGE step has to be done using the anomalous-dimension matrix up to order a z s . While 
the three-loop mixing among the four-quark operators Oi (i = 1, . . . , 6) [37] and among 
the dipole operators Or, 8 [38] are already available, the four-loop mixing of the four-quark 
into the dipole operators is still an open issue, (c) To achieve NNLL precision, the matrix 
elements (X s j\Oi(fib)\b) have to be calculated to order a 2 precision. This includes also 
bremsstrahlung corrections. In 2003, the (a 2 Np) corrections to the matrix elements of the 
operators Oi,0 2 ,Or,0 8 were calculated [39]. Complete order-a^ results are available to the 
(Or, Or) contribution to the decay width [40]. Recently, also order a 2 terms to the photon 
energy spectrum (away from the endpoint -E™ ax ) were worked out for the operator Or [41]. 

In this paper we present our calculation of the 0(a 2 ) contribution to the partonic 
b — > XP arton 7 decay rate due to the electromagnetic dipole operator Or- Contrary to [40], 
where this self-interference of the Or operator was extracted from the imaginary parts of the 
corresponding 3-loop diagrams via the optical theorem, we calculate the individual cut con- 
tributions directly. The advantage of our approach is that it also allows for the calculation 
of the NNLL interference term of the electromagnetic operator Or and the chromomag- 
netic dipole operator Og and also of the NNLL self-interference term of the operator Og 
in a straightforward manner. Applying the optical theorem in a manner analogous to the 
(Or, (^-interference is not possible in the latter cases, because some cuts contribute to 
processes other than the b —>■ XP artonic 7 decay, and they cannot be separated. We already 
anticipate here that our results are in complete agreement with the ones presented in [40]. 

The remainder of this paper is organized as follows. In section 2 we first present our final 
results, while in section 3 we discuss in some detail the various steps and ingredients of our 
calculation. Section 4 is reserved for our conclusions. 

2 Results 

Within the low-energy effective theory the partonic b — > X s ^f decay rate can be written as 
T(b^xr° n l) Ej >E = G >™<^< \ Vth V t :\ 2 ^CfMCfMG^VL), (2.1) 

where mj and ra&(//) denote the pole and the running MS mass of the b quark, respectively, 
Cf s ([i) the effective Wilson coefficients at the low energy scale, and E$ the energy cut in the 
photon spectrum. In the present paper, we calculate the function Grr(E ,{i) for E = 0, i.e. 
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without a lower cut-off for the photon energy. The function G 77 (E ,fi) corresponds to the 
self-interference of the electromagnetic dipole operator 

°7 = 7^ m M (sa^P R b) F, u (2.2) 

including a 2 terms as required for NNLL accuracy. When working to this precision this 
function can be separated into three different parts, 

4 

G 77 (0,//) = ^G^ n (0,/i), (2.3) 

n=2 

corresponding to the n particles in the final state, namely the b — > 57 (n = 2), b — > s^g 

(n = 3), b — > S'jgg and b — > s^qq (both n = 4, q e {u,d, c, s}) transitions. The individual 
contributions can be written as 

G 7 r 2 (o, /i) = 1 + ^ x 21 + 2 X 22 + 0{*l) 

G 7 r 4 (o,/i)= (^r) 2 ^ 2 +o(«D, (2.4) 

where a s (/j) is the running coupling constant in the MS scheme. The functions X{j themselves 
can be naturally divided into two parts, 

Xa = + X% , (2.5) 

where the superscripts 'bare' and 'ct' denote the contributions stemming from the unrenor- 
malized amplitudes and from amplitudes involving at least one counterterm Z factor, re- 
spectively. Working in the Feynman gauge for the gluon propagator and defining 

L, = In (2.6) 

the contributions of order a s read (inserting the color factors in numerical form) 
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while those of order a 2 are given by 
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X£ = 0. (2.8) 

These results are valid for Np quark flavors out of which Nl are massless, and Nh = Nf — Nl 
have masses equal to m&, i.e. we assume that all quark masses except mj can be neglected. 
We will address this issue once more at the end of this section. 

We note that we have numerically calculated the bare terms of order a 2 with Maple [42], 
requiring ten digits precision. However, when writing the corresponding results in (|2.8|) . we 
retained only 6 significant digits. 2 . 

When summing all up, the poles in e cancel and we obtain the final result 

G 77 (0,^) = 1 + ^(^-^-1^ 
4vr 1 9 9 3 " 
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(3.55556iV> - 44.4446)L 2 + (21.6168iV F - 334.803)L A 



+ 37.8172iV L - 2.16077iV H - 519.250 j , (2.9) 

noting that all the 6 significant digits given in this equation are in agreement with the 
numerical evaluation of the analytical result for Gtt(0, fi) given in ref. [40]. 

We also note that the terms of order a 2 Nl were calculated some time ago in ref. [39] , 
using a different regularization scheme for infrared and collinear singularities. Taking into 
account that in [39] the running mass m b (fi) appearing in (J2.1|) was expressed in terms of 
the the pole mass nib, we find perfect agreement. 

Blokland et al. [40] already investigated in detail the size of the 0(a 2 s ) corrections to 
the ratio R = T(b — > X s 7)/T(& — > X u eu) and the goodness of the naive non-abelianization 
hypothesis in various conventions for the factor of m& that normalize the decay rate. As we 
completely agree with their findings, we do not discuss this issue again. 

Instead, we briefly investigate the sensitivity of the order a 2 corrections (of course only 
those related to the (Or, Oj) interference) on the charm quark mass. As discussed above, 
the results in (12.9)1 are valid for Njj flavors with mass mj and for Nl massless flavors. As 
the charm quark mass is somewhere in between, it is interesting to consider the coefficient 
of (a s /(4vr)) 2 in (jUD for the two extreme cases (N H = 1,N L = 4) and (N H = 2, N L = 3). 
Denoting this coefficient by X2(Nh, Nl), we obtain (for \x = nib) 

X 2 (l, 4) = -370.142 X 2 {2, 3) = -410.120. (2.10) 

2 Actually, we also obtained the results for the bare terms proportional to a 2 s Nl in analytical form. 
However, to get a uniform presentation of the results, we give these terms in numerical form. 
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The true answer is expected to be somewhere in between. As a consequence the theoretical 
uncertainty of this coefficient due to the uncalculated m c -dependence is about ±5%. 



3 Details of the calculation 

As already discussed in the introduction, we do not calculate the self-interference of the 
operator O-j by taking the imaginary parts of the corresponding 3-loop diagrams via the 
optical theorem, as done in [40], but we directly work out the individual two-, three- and 
four-particle cut-contributions as indicated in (J2.4J) . 

In a first subsection we briefly summarize the main practical tools used for the derivation 
of the unrenormalized contributions to these cuts (denoted by X^ are in section^. In a second 
subsection we illustrate that we can obtain to very high numerical precision the results for 
specific diagrams in [40] when summing the contributions of the corresponding cuts. The 
relevant renormalization constants for the various fields and parameters which are needed to 
calculate the counterterm contributions can be found in the third subsection. 



3.1 Regularized bare contributions 

The starting point is the general expression for the decay rate of the massive b quark with 
momentum into 2 < n < 4 massless final-state particles with momenta fcj, 



r - = s; (n/ j^k) (»-E* ) 1*1 




\M n \ 2 , (3.1) 

-i=i J 

where the squared Feynman amplitude \M n \ 2 is always understood to be summed over final 
spin-, polarization- and color states, and averaged over the spin directions and colors of the 
decaying b quark. For n — 4 it also includes a possible factor of 1/2 if there are two identical 
particles in the final state. Furthermore, d = 4 — 2e denotes the space-time dimension that 
we use to regulate the ultraviolet, infrared and collinear singularities. 

At this point a remark concerning the polarization sums over transverse gluons is in order. 
There are two well known procedures: a) One takes the full expression for the transversal 
polarization sum, i.e 



where the polarization vectors e r , the momentum k of the gluon, and the arbitrary vector t 
are chosen such that 

e r ■ k — 0, e r ■ e r > = —5 rr >, e r ■ t — 0, r, r' — 1, 2 . (3.3) 
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Figure 1: Sample contribution to ri_> n for n = 2 (left), n = 3 (center) and n = 4 (right) 
at 0(a 2 ). The vertical line separates the original Feyman diagrams which enter the squared 
amplitude |M n | 2 . Thick lines denote the massive b quark, thin lines the massless s quark, wavy 
lines the photon, and curly lines gluons. 



b) Alternatively, one replaces the polarization sum by — and includes also (unphysical) 
ghost particles in the final state. We have checked that the two methods give the same 
results. 

The next step concerns the algebraic reduction of the quantities T\^ n being composed 
of several contributions one of which is depicted in Fig. [I] for the cases n = 2, 3 and 4 at 
0(a 2 ). Following [43,44], we replace the delta functions in the phase-space representation 
(|3.1|) by the differences of propagators, e.g., 

2m \q' s — lO q 2 + lO J 

in order to convert the phase space integrations into loop integrations (which can be combined 
with possible loop integrations already present in |M n | 2 ). In this way the systematic Laporta 
algorithm [45] based on the integration-by-part technique first proposed in [46,47] can be 
applied to reduce scalar products in the numerators and the number of propagators 3 . At 
0{af) this means that we have to reduce 3- loop diagrams, see Fig.QJ for which we use the AIR 
implementation of the Laporta algorithm [48] and a self-written one in Mathematica [49]. 

After the reduction process it usually happens that in some terms the propagators in- 
troduced via (J3.4j) appear with zero or negative power. In this case the ±z0 prescription 
becomes irrelevant and as a consequence these terms cancel. In the remaining terms we 
convert the propagators introduced via ()3.4j) back to delta-functions. 

The final task is thus to perform possible loop integrations contained in \M n \ 2 together 
with the phase space integrations. For n = 2 the phase space integration is trivial, 

ri^ = ^^ffl y (4vr) 1 -fmr|M 2 | 2 (3.5) 

and the loop integrations at 0(a 2 ) can be done numerically by means of the sector decom- 
position method [50]. The remaining loop integrations at 0(a s ) contained in \M^\ 2 can be 

3 Additional kincmatical constraints on the phase space integrations have to be implemented via additional 
delta- functions in order to be compatible with the integration-by-parts equations [44]. For example the 
implementation of a lower cut-off Eq in the photon energy spectrum can be achieved by inserting the term 

Je„ \ rn b J 

on the right hand side of l|3.1|l . 
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performed via introduction of Feynman parameters as usual. The integration of the latter 
will also be done numerically by means of the sector decomposition method when combined 
with an appropriate parametrization of the three-particle phase space, see below. For the 
treatment of the n = 4 contributions we follow [51] as far as the phase space representations 
are concerned (see also [52]). However, as the structure of the matrix elements in our ap- 
plication is slightly different than in [51], certain substitutions of the integration variables 
are necessary as described below. The (numerical) evaluation of the phase space integrals is 
then again performed by sector decomposition techniques. 

To make the paper self-contained as far as parametrizations of 3- and 4-particle phase 
spaces are concerned, we look at them in turn. To this end, we define the dimensionless 
Mandelstam variables 

11 2 

Uij...k = — o s ij...k = — + k ,• + •■• + fcfc) • (3-6) 

mr h mr h 
For the three-particle process b — > s'-yg we obtain 

Fl ^ 3 = I r(I-2) m ' d " ? / dx 1 dx 2 (l - x 1 )[x 1 (l - x 1 ) 2 x 2 (l - x 2 )]^\M 3 \ 2 , (3.7) 

where the variables y^ appearing in |M 3 | 2 are understood to be expressed in terms of the 
two integration variables, 

2/12 = (1 - Xi)(l - X 2 ), 2/23 = (1 - Xl) %2, 2/13 = ^1 • (3.8) 

For the four-particle processes b — > s'jgg and b — > s'yqq we take as our starting point (see [51]) 
r i-4 = (47r)-^mg d - 9 J dx 1 dx 2 dx 3 dx 4 dx 5 [x 1 {l-x 1 )] d - 3 



d-4 



xx 2 2 (I - x 2 ) d -*[x 3 (l - x 3 )x 4 (l - x 4 )} — [x 5 (l-x 5 )] — |M 4 | 2 , (3.9) 

where the variables yij appearing in \M 4 \ 2 are understood to be expressed in terms of the 
five integration variables, 

2/12 = (1 - xt)(l - x 2 )(l - x 3 ) 
2/i3 = (yf 3 - 2/13)^5 + Viz 

2/14 = (1 - xi){x 2 + x 3 - x 2 x 3 ) - (yf 3 - yi 3 )x 5 - y n 
2/23 = xi(l - x 2 )x A 

2/24 = Zi(l ~X 2 )(1 -X 4 ) 
2/34 = ^1^2 



Vi3 = (! - x h 



x 3 x 4 + x 2 (l - x 3 )(l - x 4 ) ± 2^X2X3(1 - x 3 )x 4 (l - x 4 ) . (3.10) 



Of course this is not the only possible parametrization since the pure phase space (13. 1|) 
with |M 4 | 2 = 1 is symmetric with respect to interchanges of the momenta fcj. Indeed, this 
symmetry sometimes can be used to avoid the appearance of the square roots yf 3 in the 
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Figure 2: Three-loop diagram (o) of reference [40]. The dashed lines indicate the 2-, 3- and 
4-particle cuts (see text). Thick lines denote the massive b quark, thin lines the massless s quark, 
wavy lines the photon, and curly lines gluons. 



squared amplitude. We remark that the formulas given above for the three- and the four- 
particle phase spaces have the crucial property that the pure phase-space integrals completely 
factorize when written in these variables. 

The expression in (|3.9|) is a convenient parametrization for the 1 — > 4 decay width for 
the cases where the denominators in |M 4 | 2 do not depend on y 13 and y^. As in [51] one can 
always remap the momenta of the final state particles in a given diagram such that 2/14 does 
not occur. In our cases, however, structures of the form (a — j/13) appear in the denominator, 
where a is a linear combination of some of the remaining yij such that (a — y\ 3 ) > 0. This is 
in contrast to the applications considered in [51] where only yi 3 appears in the denominators. 
For this reason we generalize the expressions given in eqs. (32) and (33) of [51] to our case. 
Substituting x 5 in (|3.9J) by x 5 , where 

- _ 2/13 (Z5) - Vi 3 a ~ vfs ,01 -n 

2/13-2/13 « -2/13(^5) 

we obtain the representation 

ri_> 4 = (47r)~^m^~ 9 — _ 2 _ J dx 1 dx 2 dx 3 dxidx 5 [xi (1 - £i)] d ~ 3 

d— 4 d—A d—5 

xx 2 2 (1 - x 2 ) d ~ 3 [x 3 (l - x 3 )x 4 (l - x 4 )]~[x 5 (l - x 5 )] — (a - 2/13(^5)) 
x [(a-yt 3 )(a- y ^)}^ [(1 - x s )fc 3 - y+) + (a - ^3)]^ l^ 4 | 2 , (3.12) 
in which the factors of (a — 2/13) that appear in the denominators of |M 4 | 2 get cancelled. 



3.2 Illustrative example 

As mentioned, Blokland et al. [40] calculated the order a 2 s contributions to the decay width 
by working out (according to the optical theorem) the absorptive part of the 6-quark self- 
energy. At this order the (bare) self-energy is given in terms of the 3-loop diagrams depicted 
in Fig. 1 of [40]. As an example we consider their diagram (o) and redraw it in Fig. |21 The 
diagram has a two-, a three- and a four-particle cut as indicated by the straight lines. These 
cuts represent interferences which are contained in the quantities T\^ n {n = 2,3,4) defined 
above. We denote these interferences by r^ n , where the superscript stands for diagram (o). 
As we explained in detail in the previous subsection how to calculate such interferences, we 
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immediately give the result (putting \i = m;, and working in the Feynman gauge): 
(o) _ r f°%\ n (n C A (v(») , v (o) , v (o) 



n°i n = T ( g j C F {C F - ^ J ^' + ^ + Yf>) , (3.13) 

where To denotes the decay width at lowest order in QCD. The quantities Y^, Kj , y}°\ 
corresponding to the 2- 3- and 4-particle cut, read: 

F 2 (o) = 0.333333/e 4 + 1.66667/e 3 + 2.08693/e 2 - 6.64830/e - 67.1138 

Y 3 (o) = -1.00000/e 4 - 5.00000/e 3 - 0.970923/e 2 + 72.0762/e + 410.481 

Y} o) = 0.666667/e 4 + 3.33333/e 3 - 1.11601/e 2 - 61.4279/e - 316.587. (3.14) 

For the sum = F 2 (o) + F 3 (o) + y} o) we obtain 

y(°) = 4.00000/e + 26.7798 , (3.15) 
which is in perfect agreement with the analytic result 



y(°) =4/e + 



545 5(K, . 28 2l n 1 A 2 14 4 
C(3) + — 7T 2 In 2 - 14vr 2 + — vr 4 



4.00000/e + 26.7798 (3.16) 



6 3 3 v ' 3 45 
of Blokland et al. [40]. 

3.3 Renormalization constants 

In this subsection, we collect the relevant renormalization constants necessary to render all 
bare amplitudes ultraviolet finite and fix the adopted renormalization schemes. 

First, the strong coupling constant is renormalized using the MS scheme: 

a h r = ^{£)zT S ^), (3.17) 
with 

m=1 _33_2^a^) 
3 e An 

Also the b quark mass contained in the operator O-j and the Wilson coefficient Cj are 
renormalized in the MS scheme in our calculation. In the present application we only need 
the product of these two renormalization constants for which we find 

7 ms 7 ms 1 , 4 «.(/*) / 58-4iV F 543 - 26A+A (<ti\\ n( z, ( ., q) 

Zm ^ 77 ~ 1 + si^r " {-j? 27i — ) \-^r ) + 0{as) ' (3 9) 

in full agreement with [53,54]. 

All the remaining fields and parameters are renormalized in the on-shell scheme. The 
on-shell renormalization constant for the b quark mass is given by 

*« =1 _| rWe ^(£)^ +0(< a. (3 ,o, 

ll 



For those of the gluon-field and the s quark-field we find 




47T 



) 



2 



+ 0(al) . 



(3.21) 



Finally, the field renormalization constant of the b quark in the on-shell scheme can be 
found in the literature [55]: 



4 Conclusions 

We have presented a calculation of the 0(a 2 s ) contribution associated with the electromag- 
netic dipole operator Oj to the decay width for B — > X s '-f. In contrast to a previous 
work [40], the calculational method we used is straightforwardly extendable to the calcula- 
tion of the interference of the electromagnetic and chromodynamical operator and also to 
the self-interference of the latter. Our final results are in complete agreement with those 
in [40]. Additionally, we estimated that the m c -dependence of the 0(a 2 ) contribution of the 
{Oii (^-interference, which is not worked out in both papers, leads to an uncertainty of 
about 5% (of this 0(a 2 ) contribution). 

The NNLL piece calculated in the present paper is an important ingredient of the com- 
plete NNLL prediction for the decay rate of the process B —>■ Xs'-f. According to a previous 
study, it is expected that the theoretical error of the forthcoming NNLL result will be re- 
duced by factor 2 [35] . This increases the high sensitivity of the decay B — > X^ to possible 
new degrees of freedom even further. 
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